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Surface enhanced Raman spectra from molecules (bipyridyl ethylene) adsorbed on gold dumbells 
are observed to become increasingly asymmetric (Fano-like) at higher incident light intensity. The 
electronic temperature (inferred from the anti-Stokes (AS) electronic Raman signal increases at 
the same time while no vibrational AS scattering is seen. These observations are analyzed by 
assnming that the molecule-metal coupling contains an intensity dependent contribution (resulting 
from light-induced charge transfer transitions as well as renormalization of the molecule metal 
tunneling barrier). We find that interference between vibrational and electronic inelastic scattering 
routes is possible in the presence of strong enough electron-vibrational coupling and can in principle 
lead to the observed Fano-like feature in the Raman scattering profile. However the best fit to the 
observed results, including the dependence on incident light intensity and the associated thermal 
response is obtained from a model that disregards this coupling and accounts for the structure of 
the continuous electronic component of the Raman scattering signal. The temperatures inferred 
from the Raman signal are argued to be only of qualitative value. 

PACS numbers: 85.65.+h 73.23.-b 78.20.Jq 78.67.-n 


Surface enhanced optical response of molecules ad¬ 
sorbed on metal surfaces reflects the effect of strong 
local fields created by surface plasmons[TH3] as well as 
charge transfer between the molecule and metal. mz] 
Surface enhanced Raman scattering (SERS) has become 
an important diagnostic tool for molecules at metallic in¬ 
terfaces including molecular conduction junctions. iHIS] 
(See also Refs. IT6H251 for related theoretical work). At 
metal interfaces, vibrational Raman scattering is accom¬ 
panied by a continuous background resulting from in¬ 
elastic contributions of electron-hole (e-h) excitations in 
the metal. [niiiiMiH] Both components of the inelastic 
scattering signal were recently used to determine the bias 
induced heating in a molecular conduction junction. [18] 
While it is generally recognized that the electronic 
background is affected by the molecule-substrate interac¬ 
tion, these two components of the SERS signal are usu¬ 


ally treated separately. Surprizingly, we observe (Fig. SI) 
an apparent Fano-like feature that may indicate interfer¬ 
ence between these two scattering channels. Below we 
discuss the origin of this observation. 

The measurements are carried out on single silica en¬ 
capsulated gold dumbbells (Fig. SI right inset). The 


nanosphere diameter is 95 ± 5 nm and the intersphere 
spacing prior to irradiation is ^ 1 nm. On these 
dumbbells, the quadrupolar and the binding dipolar 
plasmon resonances occur near 560 nm and 780 nm, 
respectively. [55] As molecular reporter, bipyridyl ethy¬ 
lene (BPE) is adsorbed on the gold spheres prior to en¬ 
capsulation. The dumbbells are dispersed on a silicon 
nitride membrane (20 nm thick) of the TEM grid by 
drop casting them in a dilute solution. Locations and 


geometries of the nanostructures are mapped using a 
scanning electron microscope (SEM), then Raman scat¬ 
tering measurements are done under an optical micro¬ 
scope, in the backscattering geometry, using an objective 
with a numerical aperture of NA = 0.625. The exci¬ 
tation source is a continuous wave diode laser, operat¬ 
ing at A = 532 nm, coincident with the anti-bonding 
quadrupolar surface plasmon resonance. [29] The molecu¬ 
lar vibrational Raman spectra appear over a background 


continuum (Fig. SI d), which is also seen on bare dumb¬ 
bells (Fig.[Sl^). As discussed before. [50] at low intensity, 
10 /rW/p,m^, the vibrational spectrum matches that of 
the isolated BPE molecule. Two observations are most 
significant at higher incident intensity: (a) The molecu¬ 
lar lines broaden asymmetrically to eventually coalesce 
into the single asymmetric profile, similar to Fano-type 


lineshape (Figs. SI: and d). (b) The electronic and vi¬ 
brational temperatures, and Ty, inferred from the AS 
branch of the electronic Raman scattering and S/AS in¬ 
tensity ratio of the vibrational Raman line appear to 
be different. At the highest incident light intensity Tg 
reached 580 K, while no molecular AS scattering is seen, 
which would set the apparent vibrational temperature of 
the molecule, T„, to less than 300 K. More experimental 
details and results are presented in the SI. 

The central question raised by these observations is 
whether the lineshapes shown in Figs. jSlj :, d are in¬ 
deed interference features as their Fano-like lineshapes 
suggest. Specifically: can there be interference between 
the electronic and vibrational Raman scattering path¬ 
ways that coexist in this system? The intuitive answer is 
negative, because the two pathways lead to different final 
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FIG. 1: Raman spectrum recorded on (a) a bare dumb¬ 
bell, without reporter molecules at low irradiation intensity 
(~ 10 jj.W/uni^)- (b) with reporter molecules at the same low 
intensity; (c) and (d) with reporter molecules at irradiation 
intensities ~ 50 /iW//im^ and ~ 150 respectively. 

Demarcations on (a) and (b) indicate a magnification of the 
anti-Stokes (AS) region (7x and 4.5x, respectively). The green 
traces overlapping the (c) and (d) lines are hts to the Fano 
lineshape (see SI). The right inset shows a TEM image of a 
typical dumbbell. The left inset shows the temperature ex¬ 
tracted from the fit of the AS branch of the electronic Raman 
spectrum (continuum) to the Fermi Dirac distribution. No 
vibrational AS signal is detected. 


states. The model calculations presented below yield two 
main results: (a) The aforementioned interference can ex¬ 
ist and can in principle lead to the observed lineshapes. 
(b) However, comparing the detailed calculations with 
the available experimental observations suggests that an¬ 
other mechanism, asymmetric electronic sidebands dress¬ 
ing the molecular Raman spectrum, may be dominant. 

We consider a molecule chemisorbed on a metal sur¬ 
face and subjected to an external radiation held. In the 
experimental system the metal is a nanostructure (NS) 
that promotes plasmonic enhancement, but this enters 
our discussion only implicitly, through the magnitude of 
the local electromagnetic held. The model Hamiltonian 
comprises the molecule, metal, radiation held and their 
couplings, with the electron-photon and e-v interaction 
terms treated as perturbations: 


H — Ho + Vrad + Ve-v (1) 



FIG. 2: (Color online) A sketch of the processes under¬ 

lying the Stokes (S) component of Raman scattering from a 
molecule coupled to a metal, (a) The incident photons induce 
the virtual transition displayed in red and a corresponding 
Rayleigh scattering (the red process in reverse). In addition 
two inelastic scattering (Raman) pathways are shown: The 
green process brings the molecule back to its ground state 
while creating an e-h excitation in the metal (electronic Ra¬ 
man scattering). The blue process brings the molecule to an 
excited vibrational state on the electronic ground state (vi¬ 
brational Raman scattering). These Raman scattering path¬ 
ways are coupled by the combined effect of the metal-molecule 
electron- transfer coupling and the electronic-vibrational (e- 
v) interaction (orange) that induces vibrational de- excitation 
in the molecule with e-h formation in the metal. The dashed 
lines connect coupled processes, (b) Two pathways are de¬ 
fined by the excitation process (red): A (virtual) molecular 
transition and a molecule-to-metal charge transfer transition, 
respectively denoted M and G in the text. The outgoing pho¬ 
ton (blue) leaves the molecule in the excited vibrational state, 
possibly with e-h excitation (or de- excitation) in the metal 
(black). Note that the processes shown are accompanied also 
by electronic Raman scattering from the bare metal. 


where (fi. = fcs = e = 1) 


Ho — Hm + Hj^s + ^NS + H^ad] ( 2 ) 

Hm — ^ ^ ~b UJyV^V^ (3) 

m—g,x 

^NS ^ ^ ; ( 4 ) 

k 

Vns = E E( VmkSla^k + Vkmc\,dm'j ] (5) 

m=g,x k 

Ve-v= X! ^mQvd'^md7n]Qv=V + v‘^ (6) 

m—g,x 

Hrad — ^ ^ (^) 

a6i,{/} 

Vrad = Vr^ad + (») 


J=C,M aei,{f} 


Hm is the molecular Hamiltonian, modeled as a two elec¬ 
tronic levels (ground and excited states, g and x) system 
with a single molecular vibration of frequency ujy. Ve-v 
is the corresponding polaronic-type e-v coupling. Hns 
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describes a free electron metal and Vns is the molecule- 
metal electron transfer coupling, (dm) {'m=g^ x) and 
Cfc (cfc) create (annihilate) electrons in the molecular or¬ 
bitals and the metal, respectively. (D) anddi, (da) 
create (annihilate) vibrational excitations and photons 
in mode a, respectively. Finally, Dm = d^d^, and Dc = 
d^gCk are respectively intramolecular and charge-transfer 
electronic de-excitation operators. Vrad, Eq- (|^, is the 
system-radiation field coupling in which the coupling ma¬ 
trix elements connect the molecular ground state 

g to excited states x in the channels J = C or M that 
determines the nature of the excited states x (a simi¬ 
lar term, not shown, is associated with light scattering 
from the bare metal). As usual we distinguish between 
the incoming (occupied) photon mode i and the contin¬ 
uum of accepting (empty) modes {/}. These coupling 
elements are assumed to be functions of the vibrational 
coordinate, described here by their Taylor expansion up 
to linear term (J = M, C) 

u^.JQv) ^ c/i,;.(4)«+ U^lo^Qv (9) 

To account for the observed dependence of the inelastic 
spectrum on the incident light intensity, we assume that 
in addition to being scattered, the incident radiation field 
affects the molecule-metal coupling. This can arise from 
a field-induced renormalization of the electron transfer 
coupling Eats, Eq. ([^. [5TH55] as well as contributions 
from molecule to metal charge transfer transitions. In 
the calculations reported below this is incorporated by an 
assumed dependence of the self-energies T^ (rn = g, x) 
on the incident light intensity. 

Consider first model A (Fig. |S2^ ). Here we assume 
that Raman scattering is dominated by the intramolec¬ 
ular optical transition and disregard the corresponding 
contribution of the optical charge transfer coupling 
This model can show interference between the vibrational 
and electronic scattering channels, provided that their 
coupling, Eq §, is strong enough. 

Following the procedure described in Ref. [TS] the (nor¬ 
mal) Raman scattering flux has the following struc¬ 
ture 

/ +00 pO pO 

d(t' -t) d{ti -t) d{t2 - t') 

-oo ^ —oo J —oo 

X ( 10 ) 

X (u,{t2)D{h) Uf{t')b\t') Uf{t)b{t) 

in which the operators are summed over all relevant con¬ 
tributions from Eqs. Q and ^ and where the indices 
i and / indicate that the corresponding matrix element 
should be taken with the incident or scattered radiation 
field modes, respectively. In departure from Ref. [18] we 
keep only the lowest order (up to 2”“^) terms in U^, as 
usually done for non-resonant Raman calculations. The 
calculation proceeds as follows (See the SI for details): 


(a) Applying Eq. (S2) to model A, yields Eqs. (S2) 
for the light scattering flux. Analysis based on energy 
conservation suggests that expressions (S2a), (S2f) and 
(S2k) contribute to Rayleigh and pure electronic Raman 
scattering, while expressions (S2g) - (S2j) contribute to 
vibrational Raman scattering (possibly dressed by e-h 
excitations). dU Expressions (S2b) - (S2e) formally rep¬ 
resent interference between vibrational and electronic- 
Raman/Rayleigh scattering. At the level of the present 
calculation, fourth order in the coupling to radiation 
field, they vanish in the absence of the coupling([^. 

(b) The e-v interaction ®> is considered to the lowest 
order. This introduces an additional term 


J dTyMgQy{Ty)hg{Ty) (H) 

into the correlation functions (S2b)-(S2e), yielding non¬ 
vanishing corrections that represent interference between 
the electronic and vibrational Raman processes that can 
give rise to Fano lineshape in the Raman scattering as 
shown below. [45] In evaluating these terms, the real time 
analog, ty, of the contour variable Ty can be placed in all 
possible ways between the times t, t', ti, and ^2 on the 
Keldysh contour (see Fig. S5). 

(c) The light scattering flux including the interfer¬ 
ence corrections is thus obtained in terms of correla¬ 
tion functions of a system defined by the quadratic 
Hamiltonian Hq. Thus Wick’s theorem applies and 
electronic and vibrational degrees of freedom decouple 
yielding the final expressions in terms of projections to 
real time of the electronic and vibrational Green func¬ 
tions, Gmm'{T,T') = -i{Tcdyn{T) and D{t,t') = 

—i{TcQy(T) Qyir')), where is the contour ordering op¬ 
erator. For the molecular vibration we use the quasipar¬ 
ticle approximation, Eqs. (S3) and (S4). In evaluating 
the electronic GEs we assume £x — where F^ = 

|EfcrriP^(E — Efc); m = g,x. Consequently, inter¬ 
state correlations induced by the metal are neglected, 
~ GrfijYi{T^ T ) — dya^m' G X ), lead 

ing to the standard expressions (S5)-(S7) for the elec¬ 
tronic GFs. 

This procedure leads to explicit expressions, Eqs. (S8)- 
(S14), for the different contributions to the light scat¬ 
tering flux in the relevant order (fourth) in the system- 
radiation field coupling including the interference correc¬ 
tions, Eqs. (SIO) and (S14), arising from the e-v cou¬ 
pling. The other contributions can be classified accord¬ 
ing to their physical origin that can be identified by their 
energy conservation structure. Eq. (S8d) and (SlOd) 
represent the the Rayleigh scattering component that 
is disregarded in our calculation. Eq. (S8) represents 
the pure electronic Raman that may be though of as an 
electronic sideband of the Rayleigh peak.|46| We have 
previously [Ml IM] used the Stokes (S) and anti-Stokes 
(AS) components of this contribution to estimate the 
apparent electronic heating in non-equilibrium molecu- 
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FIG. 3: (Color online) (a) The total Raman scattering 

(fnll line; blue) evaluated from Eqs. (S8)-(S14) (without the 
Rayleigh contributions, Eqs. (S8d) and (SlOd)), displayed 
against the Raman shift Vi — Vf. The parameters nsed are 
T = 300 K, eg = 0, = 2 eV, Tg = T, = 0.05 + eV 

with = 0.3 eV being the light induced tunneling rate, 

= 0.05 eV, Mg = 0.02 eV, Vi = 1 eV, = 0.1 eV, = 
0.01 eV. The red dashed line is the resnlt obtained when 
the e-v interaction ([^ is disregarded. The inset shows the 
corresponding weak incident held case with = 0.03 eV, 
= 0.032 eV, = 0.0032 eV. (b) The Raman scat¬ 
tering calculated with the same parameters as in the main 
panel on the left at temperatnres 3000 K (solid and dashed 
lines) and 30 K (dotted and dash-dotted lines, mnltiplied by 
a scale-factor 4). 


lar junctions. Eq. (S13) is the vibrational contribution 
to the Raman signal together with its own electronic side¬ 
band. 

Fig. shows the total Raman scattering signal re¬ 
sulting from this model calculation. The calculation was 
done on an energy grid spanning the region from —3 to 3 
eV with step 10“^ eV and the molecular vibration is as¬ 
sumed to equilibrate quickly to the ambient temperature. 
The following observations can be made: 

(a) Interference between electronic and vibrational in¬ 
elastic scattering routes can indeed lead to a Fano-type 
scattering spectrum. The electron-vibration coupling ®> 



FIG. 4: The effective temperature obtained from the Raman 
scattering spectrum displayed in the main panel of Fig. 3. 
See text for details. 


is crucial for obtaining this behavior. This is seen in the 
weak illumination case (inset to Fig. S3 r) by comparing 
the results obtained with and without this coupling. 

(b) The asymmetric lineshape that seen in the high 
incident intensity case (main hgure S3i) is not of this 
origin. In this case the calculation yields Fano-like line- 
shapes that qualitatively reproduce the experimental be¬ 
havior also when the coupling (|^ is disregarded. 

(c) The transition from a standard, essentially sym¬ 
metric Raman line to a broadened assymetric line when 
the incident light intensity increases is best rationalized 
within this model by assuming that the coupling is 
small enough to be disregarded. This yields the dashed 
lines in Figure that qualitatively agree with the ex¬ 
perimental observation. 

Our analysis thus shows that interference between vi¬ 
brational and electronic Raman scattering can in prin¬ 
ciple occur leading to a characteristic Fano lineshape in 
Raman scattering from molecules adsorbed on metal sub¬ 
strates. It suggests, however, that the origin of of the 
present observation of asymmetric Raman lineshape at 
higher (photoinduced) metal-molecule coupling is differ¬ 
ent resulting from the S component of the electronic Ra¬ 
man scattering. This is seen in the Fig. |S3^ 3, that show 
that the asymmetry disappears at high ambient temper¬ 
atures. 

Our calculation also yields the effective (Raman) tem¬ 
perature from the ratio between the S and AS scat¬ 
tering intensities according to [131 dZl [T5] = 


Ai/ 


ln( 




, or htting the tail of the AS sig- 

(2)i 


nal to a IS.v /(l — exp[—Eq.(Sl), where IS.v = 


\vi — i'f \ is the Raman shift. For we have arbi¬ 

trarily assigned the vibrational Raman temperature to 
the value calculated at = oj^ and the electronic Ra- 
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FIG. 5: The Raman scattering signal calculated for model B 
(see text for parameters). The dotted-black and dashed-red 
lines represent the signals associated with the pure-electronic 
and electronic-dressed vibrational Raman, respectively, and 
the full blue line is their sum. The inset shows the corre¬ 
sponding vibrational and electronic effective temperatures, 
calculated as before. 


man temperature to that obtained at lS.v = 2uiy. The 
results for the parameters of Fig. S3 1 (main panel) are 
shown in Fig. |S4[ The resulting trends are similar to 
those observed experimentally, although the model cal¬ 
culation somewhat underestimates the electronic heating. 
In spite of this qualitative agreement with the experimen¬ 
tal observation, we note that while can give a rough 
estimate of system heating, the actual numbers may be 
meaningless. This is indicated here by discrepancy be¬ 
tween the different calculations of T^i and the apparent 
cooling seen in Tyn,, keeping in mind that the oscillator 
was kept at the ambient temperature in this calculation. 

We have also considered model B as an alternative 
mechanism for interference. The calculation follows the 
same steps as in (a) and (c) above, now using both M and 
C terms in Eq. ([^ and again keeping only terms up to 
the lowest order in (Eq. ([^). In this calculation 
we have restricted our considerations to contributions 
from the purely intra-molecular processes {D = Dm in 
all four terms in (S2)), purely charge-transfer processes 
(D = Dc in all terms), and their mixtures (potential 
source of interference) in which excitation operators at 
times t 2 and t' are of one kind, while those at times t 
and ti of the other. Again the molecule-metal tunnel¬ 
ing coupling Fm is assumed to contain an incident field 
dependent contribution. m The final result for the light 
scattering flux for this model is given by Eqs. (S15)- 
(S19). 


Fig. S5 shows the Stokes side of the Raman response 
obtained from this model using the following parame¬ 
ters: T = 300 K, Eg = 0 eV, Sx = 2, eV, = 0.1 eV, 


Vi = 1 eV, Fg = F,, = 0.15 eV-br°P‘ with = 0.1 eV, 

U^Ma — 0-1 eV and U^Ma — “ b/)j 

Ep = 0. Calculations are performed on an energy grid 
spanning the region from —5 to 5 eV with step 10“^ eV. 
Asymmetric Raman profiles are seen as before, but again 
we find (see SI) that while there is some contribution 
of interference between channels M and C, the main 
source of lineshape asymmetry originates again from the 
electronic sideband that dresses the Rayleigh and vibra¬ 
tional scattering peaks mainly on their blue side. The 
“electronic heating” obtained here is stronger than in 
model A, closer to the magnitude observed experimen¬ 
tally, but the actual numbers should again be considered 
cautiously. 


In conclusion, we have observed asymmetric lineshape 
features in Raman scattering from bipyridyl ethylene 
molecules adsorbed on gold nanostructures, with the fol¬ 
lowing characteristics: (a) Asymmetry increases with 
incident light intensity, (b) The vibrational tempera¬ 
ture appears not to increase even at the highest inten¬ 
sity used, while the apparent electronic temperature in¬ 
creases by up to ~ 600 K. The Fano-like appearance of 
these lineshapes suggests the possibly implication of in¬ 
terference between different scattering pathways. Model 
calculations show that such interference, leading to the 
observed asymmetry, is indeed possible in systems with 
strong enough electron-vibrational coupling. Our cal¬ 
culations suggest, however, that the observed lineshape 
asymmetry is dominated by electronic scattering side¬ 
bands that dress the Rayleigh and vibrational scattering 
peaks, and its dependence on incident light intensity can 
be explained by an optically induced component in the 
molecule-metal electron-transfer coupling. This model 
also yields vibrational and electronic “Raman tempera¬ 
tures” that are consistent with experimental observation, 
but may reflect the complex nature of the non-equlibrium 
response rather than the actual temperature. Fano-type 
interference in Raman scattering has been shown to be a 
theoretical possibility. 
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SUPPLEMENTARY INFORMATION 


Experimental Details 


s outlined in the main text, the measurements are carried out on single silica encapsulated gold dumbbells, the 


transmission electron micrograph (TEM) of which is shown in the inset to Fig. SI The nanosphere diameter is 95±5 
nm and the intersphere spacing prior to irradiation is ~ 1 nm. As molecular reporter, bipyridyl ethylene (BPE) is 
adsorbed on the gold spheres prior to encapsulation. The dumbbells are dispersed on a silicon nitride membrane (20 
nm thick) of the TEM grid by drop casting them in a dilute solution. After mapping out locations and geometries 
of the nanostructures using a scanning electron microscope (SEM), Raman scattering measurements are carried out 
under an optical microscope, in the backscattering geometry, using an NA = 0.625 objective. The excitation source is 
a continuous wave diode laser, operating at A = 532 nm, which is resonant with the anti-bonding quadrupolar plasmon 
on these structures. |29L [ 86 ] On these dumbbells, the binding quadrupolar plasmon and the binding dipolar plasmon 
resonances occur near 560 nm and 780 nm, respectively. |29j The molecular vibrational Raman spectra invariably 
appear over a background continuum (Fig. a, which can also be seen on bare dumbbells (Fig. la). The molecular 
lines broaden asymmetrically as a function of irradiation intensity, to eventually coalesce into the asymmetric Fano- 
like profile shown in Fig. |S1| ; and Id. We see variations in details on different particles, along with hysteretic response 
during intensity cycling due to evolution in the structure of the intersphere junction, which is verified through TEM. 
A rich variety of phenomena can be identified during fusion of the plasmonic junction when the incident light intensity 
increases. [3S] Here, we focus on the development of the Fano-like profile that possibly indicates interference between 
the electronic Raman scattering continuum on gold and the discrete resonances of the molecule. The temperature of 
the molecular vibrations is inferred from the usual ratio of Stokes/Anti-Stokes scattering. The electronic temperature 
is obtained from the anti-Stokes branch of the continuum, which for Raman shifts £ = oJi — iOg > kT, can be seen 


in Fig. S2 to decay exponentially, e:xjp{—£/kTf,), as expected for e-h Raman scattering that terminates on thermally 
populated holes. The spectrum is explicitly fitted to the joint density of states, which appears in first term of Eq. (SI I. 
The principal experimental observations are summarized as follows: 


1. At low intensity, 10^W//im^, the molecular vibrational Raman spectrum that appears over the background 
continuum in Fig. perfectly matches that of the isolated BPE molecule. Note that the traces in Fig. |Sl|are 
identical to those in Fig. 1 of the main text. This has been shown in some detail previously. [50] 


2. The background continuum, which is also present on the bare dumbbell (Fig. |Sl^ ), can be assigned to electronic 
Raman scattering (ERS) on the gold nano antenna. The continuum is polarized. As in the Hertzian antenna. |87j 
the scattering is dipolar, with cos^q polarization along the long axis of the dumbbell. 


3. The fits of the anti-Stokes spectra to Eq. SI I, which are shown in Fig. S2 yield the effective electronic tempera¬ 
ture of the gold. This has been independently recognized recently.^ Note, the method does not involve a ratio 
between Stokes and anti-Stokes scattering, which can be misleading. Te; is a linear function of the irradiation 
intensity, as seen in Fig. I 


4. As the excitation intensity is increased up to 400 /rW//xm^, the molecular lines develop asymmetric profiles 
and gain intensity relative to the background. At high intensity, the molecular response collapses into a single 
asymmetric line. The spectral profiles [ST}; and |Sl[ l can be fitted to a sum of ERS background and a Fano line, 
according to 


W{£)=a J f{E) (^1-f{E + £)'^e-^^^/‘^^dE + b 


1 + e'^ 


- 1 


(SI) 


in which a,b are normalization constants (a/b =20), f(E) is the Fermi-Dirac distribution, d = 0.19 eV is the 
energy width of the projection of the collective plasmon state on single particle states, [33] e = is the 

reduced frequency with center on the molecular vibrational frequency, 'Wo= 1600 cm“^. The resulting fits are 
shown as green traces overlapping the lines c and d in Fig. SI In Fig. SI: the fit parameters are 7 = 130 cm“^ 
and q = 2.5, while in Fig. Sll 7 = 390 cm“^ and q = 2.7. 


5. While the ERS background shows linear dependence on irradiation intensity, the integrated area under the 
molecular lines is superlinear, as shown in Fig. |S 3 | 3 . 
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FIG. SI: (a) Raman spectrum recorded on a bare dumbbell, where reporter molecules are absent at the metallic junction, 
(b) Molecular Raman lines of BPE appear over the background continuum from a dumbbell. Both (a) and (b) have been 
recorded at low intensity (~ 10/rW//im^) and anti-Stokes region is expanded vertically for clarity in either case (7x and 4.5x, 
respectively). At elevated irradiation intensity the entire spectrum develops Fano-like profile. Parentage of the molecular lines 
may still be identified (c) at moderate intensity (~ 50/iW//im^) and becomes unidentifiable (d) at more intense irradiation 
(~ 150/iW//rm^). The entire line profile in (c) and (d) can be fitted to Fano lineshape (green trace). The inset shows a 
TEM image of a typical dumbbell. Shaded regions are used to quantify the intensity dependence of molecular lines (blue) and 
electronic Raman scattering of gold (red), the results of which are shown in Fig. SI d. 


6 . Remarkably, the molecular vibrations and metal electrons appear to have different temperatures. At the highest 
incident light intensity, where the apparent electron temperature reaches = 580 K, we see no evidence of 
molecular anti-Stokes scattering. Based on the Stokes to anti-Stokes ratio, we establish that the apparent 
vibrational temperature of the molecule, T„, is less than 300 K. 

7. Near 1 mW/^m^, the nano-structure undergoes phase explosion - the metal melts and the silica shell explodes. 
Hysteretic behavior is observed during intensity cycling at intensities in excess of ~ 100/iW//rm^, due to per¬ 
manent evolution of the junction structure. The spectrum in Fig. |S1| :, which was recorded at an irradiation 
intensity of 50/rW/^m^, was obtained after such a cycle. 


Theoretical Details 


Here we present details on derivation of Raman signal expressions for the two models considered in the paper. 
The starting point is general expression of the normal Raman scattering (see Eq.(9) of the paper and Fig. |S4| below) 

















FIG. S2: The anti-Stokes branch of the continuum is given by the thermal occupation of the Fermi-Dirac distribution. This is 
illustrated by the exponential fits, Ias{£) = exp(—f/feTe), to the anti-Stokes continuum for four different irradiation intensities, 
I = 100, 50, 25,13/iW//im^, which yie ld e lectronic temperatures of Te — 465,415,395,364 K, respectively. A more accurate 
value of Te is obtained from fitting to (SI I, but the trend does not change. 
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Intensity (|jW/|jm^) Intensity (jjW/|jm^) 


FIG. S3: (a) Exponentially decaying anti-Stokes branch of Raman spectra determines the temperature at the metallic junction 
which grows linearly with increasing irradiation intensity, (b) The electronic Raman scattering background (red dot) varies 
linearly with irradiation intensity. However, the area under the molecular lines (blue dot) is superlinear. 
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FIG. S4: (Color online) Projection for ‘the normal Raman’ scattering. Times ti and t 2 indicate interaction with the incoming 
photon Vi, t and t' - with the outgoing photon vj. The diagram is identical to Fig. 8a of Ref. |18 |. 


[M SOI mi- 

Utilizing Taylor expansion of the optical electronic transition matrix element (see Eq.(8) of the paper) and keeping 
terms up to second order in molecular vibration coordinate leads to the following steady-state expression for Raman 
scattering from mode i to mode / of the radiation held 


/ +00 1*0 /*0 

d(t'-t) / / d{t2 - t') - iv^(ti 

-oo J —oo J —oo 

D\t') D{t) Dt(U)) 

+ D\t') b(t) bHh)) 

+ u\%U^^]jUfluflbD{t2) bbt!)Ql{t') m D\h)) 

+ bbf) mQ,{t) bbu)) 

+ bbf) b{t) bHtbQlih)) 

+ bbt')Ql{t') m Dbh)) 

+ ul%u^gl]fUj%ull]bD{t2) bbt')QUt') bit) bbh)Qlih)) 

+ u^%u^Zu^^%u^^l]bD{t2) bbt') bit)Q,it) bbtbQliti)) 


^2) 


(S2a) 

(S2b) 

(S2c) 

(S2d) 

(S2e) 

(S2f) 

(S2g) 

(S2h) 

(S2i) 

(S2j) 

(S2k) 


Analysis based on energy conservation suggests that from all the resulting diagrams expressions (S2a), (S2f| and 
( |S2k[ ) contribute to electronic Raman scattering (in particular, Eq. ( |S2a[ ) leads to results discussed in out previous 
publications [5il 155] ). Expressions (S2g)-(S2j) contribute to vibrational Raman scattering (these results are similar 
to those of resonant Raman scattering consideration in our previous publications [171 [18]). Einally, expressions (S2b)- 
( S2e[ ) are responsible for interference between the electronic and vibrational intra-molecular Raman scattering channels 
(see details below), which is possible in the presence of electron-vibration interaction (see Eq.(5) of the paper). This 
situation is encountered in description of the model A (see Fig. 2a of the paper). If one disregards this interaction (and 
within the treatment of molecular coupling to radiation held up to second order in molecular vibration coordinate) 
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the latter contributions are zero. This is the situation encountered in description of the model B (see Fig. 2b of the 
paper). We now turn to details of derivation specific for each of the models. 


Model A 


Within this model we are interested in interference between electronic Raman scattering (a process where vibrational 
state at the start an end of Raman scattering is the same) and vibrational Raman scattering (a scattering process 
with initial and final states with different vibrational excitations) amplitudes. For the two amplitudes to interfere, 
final states of the two processes should be the same. Electron-vibration coupling in molecular Hamiltonian (see Eq. (5) 
of the paper), yields such a possibility (see Fig. 2a of the paper). In the off-resonant Raman regime electron-vibration 
interaction can be disregarded. Utilizing the first order perturbation in the electron-vibration coupling (see Eq.(lO) of 


the paper) in expressions (S2b)-(S2e) results in the desired interference. Possible placements of the time of electron- 


vibration interaction between times t, t' and U, ^2 representing interaction with the outgoing / and incoming i 
modes of the radiation field, respectively, are shown in Fig. [S^ 

Projections (e)-(h) of Fig. S5 utilized in corrections to ( S2b| )-(S2c) and projections (a)-(d) of Fig. S5 in corrections to 
(S2d)-(S2el lead to renormalization of electronic Raman due to electron-vibration interaction. These are the diagrams 


which lead to Fano resonance in Raman scattering (see below). Projections (a)-(d) of Fig. in corrections to (S2b)- 
(S2cl and projections (e)-(h) of Fig. S5 in corrections to (S2d)-(S2e) renormalize vibrational Raman scattering. 


After accounting for the electron-vibration interaction we decouple electron and vibration degrees of freedom and 
utilize the Wick’s theorem to get final expressions in terms of projections of electron and vibrational Green functions 
(see Eqs. (11) and (12) of the paper). We treat molecular vibration within quasiparticle approximation 


D^{uj) = — 2TTi {Ny6{u} — (jjy) + [1 + Ny]d{uj + ojy)) 

D^(uj) = — 2TTi {NyS{uj + ujy) + [1 + Aly](5(a; — uiy) 


(53) 

(54) 


where Ny is average population of the vibration. Explicit expressions for electronic Green functions are (m = g,x) 

iTrufiE) 


G^{E) = 
G^{E) = 


(F;-e„)2 + (r^/2)2 

-tTmfiE) 


(E-e^)2 + (r^/2)2 
Gl,{E =[E -e^ + *F„/2]-i G^^^iE) = [G^^iE)]* 


(55) 

(56) 

(57) 


After tedious but straightforward algebra we get the following steady-state expressions for Raman scattering. 



FIG. S5: (Color online) Electron-vibration interaction induced corrections to ‘the normal Raman’ (Fig. |S4[ ). indicates time 
of interaction with molecular vibration. The contour is deformed following the Langreth rules [3^ . 
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Explicit expression for the electronic Raman is 


Here 


fdEgl fdEg2 fdE^l fdE,2^ ,, Z. Z. ^ 

J ^ 


G<{Eg,)G>{Eg2)G<{E,,)G>{E,2) 


+ 2 ReC/^ic/S;C/S,C/S, 


ry-(O) Ty-(O) 

^ f,gx'^gx,i 

Vi — Ex2 + Egl + i5 
l + iV„ 


(^-i - E^ 2 gi - a;„ + i(5)(r'i - Ex2gi - id) 
N„ 


{vi - E^2gl +(^v+ iS){vi - E^2gl “ *<5) 

+ J J 2TTS{l^if - Eoc2i) {Exi)G^ {Ex2) 

"dEg c/}K!zG<(e; '' 


2tt Ui — Ex 2 + Eg + id 


x 2 R*t/™Z/™,Z/™,Z/™ 


JiE^J 


1 + Ny 


{vi - E^2gl -UJy+ id){iyi - E^2g2 - id) 
N„ 


+ 


{vi — E^2gi + iOy + id){vi — E^2g2 ~ id) 
j J - Eg2i)Gf{Egi)G^{Eg2) 


dE, ufluflG>{E,) 


27r Vi — Ex H“ Eg\ -\- %6 




dExi f dEx2 

27 r J 27 r 

1 + iV, 


G>iEy,)G>iEy2) 


+ 2'!rd{h'i — Vf) 


{vi - E^igi -ujy+ id)(vi - E^ 2 gi - id) 

^ _ Nv _■ 

{vi - E^igi +ujy+ id){iyi - E^2gi - id) 

dEg f dE,uflufJ^fi<{Eg)G>{E,)^ 


27r J 27r — Ex E Eg -\- i5 


fd^ fd^ [•^3aa<{Ey)G<{E,,)Gl(E„)a>jE..,) 


27T 


27T 


271 


271 




+ 


g 

1+Ny 


{vi - E^igi -ujy + id){vi - E^2g2 - id) 
Ny 


(vi - E^igi +UJy+ id)(v^ - E^2g2 - *(5) J 


(S8a) 


(S8b) 


(S8c) 


(S8d) 


Vif — Vi Idf Exnij^gmg — Exm^: Egyy^ E^21 — Ejji2 Ej^i {tTI — x) 


(S9) 
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Corrections to the electronic Raman due to electron-vibration coupling are 




dE. 


'si 


dE. 


S2 


dE^i 

2TT 


dEx2 

2TT 


2tt J 2tt _ 

2^6{v,f - Eg2i - E,2 i) GfiEg,)G>{Eg2)G<iE^,)G>iE^2) 


X 2Re 


N„ 


AO) tAE 


r(o) tAo) 


■N„ 


— i 


{yi Ex2gl E)x2gl ^ 

iyyyijGl{E,2 T UJ.) + T 


\i^i - E, 


x2gl 


■Z,5|2 


dE, Gf{E,,)ufiyy (Uyu'^;'Jj[E{E,2,Eg,) - E{Eg2,E,,)] 


27r 


r(i) rr( 0 ) 


Vi — Ex2gl -\- i6 \ Vi — Ex2gl — Eg2g3 — iS 

r( 0 ) r^l) 


UlZyjjinEgs, E,,) - F{E,,, E,,)] 


^ PtlKJjHEgi) + UZUlxy*iEg2)]U}: X’, 

Wi - E^2gl + i5\^ 

- 2TT5{v,f - E,2 i) G<{Eg,)G<{Eg2)G<{E^,)G>{E^2) 


AO) tAA 


Ex2g3 

(0) rr(0) 


X 2Re 


1 + [Utyx,fGl{Eg2 T + UtgxU^]fG'^giEg2 T iOv)]uf!luit 

Ey J {Vi Ex2gl H“ Ex2g2 ^ 


•M 


[ul%u‘'yGl{E,2 T u:x) + UyuflfGl{E,2 T 
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{Ui - Ex2gl + iS){vi - Ex2g2 “ iS) 
rio) ttW ^WTj)o) jj{o) -, 

gx,i 


— 2TrS{h'if — Eg 2 i) G^{EgijG^{Eg 2 )G^{Exi)G^{E x 2 ) 2Re 


Uyuf]jGl{Eg2 T ^ ul%ufljG-{Eg, t 


Vi — E, 


xlgl 


-iS 


1 + iV, 

N„ 


jAo) jAo) 


■Nx 


u(o) Tjio) 

O’ f.gx'-^gx.l 


Exlg2 E J Ex2gl ^ 

uSxGixjGliEg2 T ^x) ^ UyU^^JjG’^giEg, t u^x) 


Ex2gl H“ i(5 Y Vi Exlgl ^ UJy Exlgl iS 

■ dE, Gf{Eg,)uflufl (UyufyF{Eg 2 ,Eg,) - F{Eg 2 ,Eg,)] 


2'k Vi — E. 


xlgl 


iS 


lESxG^^JyiEg,) + u^^;;^iu^^’j<i>*(Eg2)]u}';’y:’, 


+ 

AO) rr(i) 


— Ex2gl — Eg2g, — i6 

uEGgl]fiPiEg3,Eg,) - FiEg2,Eg,)] 


Vi — Ex2g3 — 
( 0 ) TtE) -\ 


{v^ - Exigi + i5)(vi - Ex2gi - id) 

+ 2ti5{v, - vf) G<iEg,)G<{Eg2)G>{Ex3)G>{Ex2) 


1 + Nx\ PSxU^x]fGl{Eg2 T uix) + U^XjGgiEg2 T uix)]U}'^>xU^Z 
N„ 


AO) rr(l) 


r( 0 ) rdO) 


■N, 


iyi -£* 3:151 E iE){yi Ex2g2 E^v 

PEu'gxjGl{Eg 2 T iOx) + uEuf]jG^g{Eg 2 T o^x)]ufluE 

{vi - Exlgl + i6){vi - Ex2g2 - iS) 

PSEtE^O^) + GEGglE*(^92)PEuE ' 

[vi - Exlgl + i5){Ui - Ex2g2 - id) 


(SlOa) 


(SlOb) 


(SlOc) 


(SlOd) 


X 2Re 
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where upper (lower) row and sign correspond to Stokes (anti-Stokes) scattering channel, repectively, and 


F{Ei,E2) = 


1 


^E) ^ z / ^ 


— £'21 

G>iE') 


+ E 21 — id 


G<{E') 


2tt \ E — E' — (jjy — iS E — E' + ujy — id 


2G<{E'] 

IjJv 


(511) 

(512) 


At steady-state explicit expression for the vibrational Raman is 


dE. 




dE, 


S2 
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27r 


dExi 
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dEx2 E 
2 tt I Nx 
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2Re 
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FIG. S6: (Color online) Intra-molecular Raman scattering (model A, see Fig. 3 of the paper). Shown are the total signal (blue 
line), electronic Raman (magenta) Eqs. (S8a|-(S8dl and its corrections due to electron-vibration interaction (red) Eqs. ( S10a| |- 
( SlOdl, as well as vibrational Raman (black) Eqs. (S13a|-(S13dl and its corrections due to electron-vibration interaction (green) 


Eqs. (|S14al-(S14dl. 
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Corrections to the vibrational Raman due to electron-vibration coupling are 
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(S14b) 


(S14c) 


(S14d) 


Similar to our previous considerations [841185j we are interested in Raman scattering only. Thus we drop Rayleigh 
contributions, Eqs. (SSdl and (SlOd). The latter are sharply peaked peaked at Vi = Vf. 

Results of numerical simulations presented in Fig. 3 of the paper show asymmetric scattering profile that becomes 
more pronounced at higher incident light intensity. Interference between the electronic and vibrational scattering 
pathway contributes to this asymmetry as seen in Fig. |S6| This interference is associated with the corrections to the 


electronic Raman due to electron-vibration interaction, Fqs. (SlOa)-(SlOd) (the orange line showing the correction 


due to the coupling M to the electronic Raman scattering in Fig. S6). However, the dominant contribution to the line 


asymmetry is seen to be the Stokes nature of the electronic Raman scattering (purple line in Fig. S 61 that dresses the 
vibrational peak on its high energy side. 


Model B 


Within the model B (see Fig. 2b of the paper) we focus on intra-molecular and charge transfer contributions, 
and disregard electron-vibration coupling (see Fq.(5) of the paper). Thus expressions (S2b)-(S2e), which are the 


source Fano-type interference in model A, do not contribute in this model. Within the model radiation field leads 
to optical excitation of either intra-molecular or charge transfer character. For each type of the excitations we 
separate vibrational and electronic degrees of freedom in the correlation function, disregard coherence between ground 
and excited states of the molecule, and perform integrals. A lengthy but straightforward derivations yield explicit 
expressions for electronic and vibrational Raman signal. The steady-state Raman flux from intial mode i to hnal 
mode / of the radiation field is 
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FIG. S7: (Color online) Vibrational (top) and electronic (bottom) Raman scattering (model B, see Fig. 5a of the paper). 
Shown are the total signals (dashed red line - top graph, dotted black line - bottom graph), as well as intra-molecular (solid 
black lines) and charge-transfer (dotted magenta lines), and their interferences (dash-dotted green lines). 
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Steady-state Raman flux from initial mode i to final mode / of the radiation field is 
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- Ex ^ ^ + ll5 G - ^ “1“ E/g ^ “1“ “1“ 


■ 2Re ^ G<(i?(1))G< {Ef))g>{ eA))G> (4^)) 


rr(l)rr(0)./l)..(0) 

Mf'v fC'^Ci 


^ {vi — Ex ^ + Eg ^ + id){vi — Ex ^ + Eg ^ + sw^, — id) 

rjAhj(O) rAO).j(l) 

_l_ Mf'v fC'-^Ci _ 

[vi — Ex ^ + Eg ^ + sWt, + id)[vi — Ex ^ + Eg ^ + sw^, — id) 

_l_ ^iM'-^ Mf^ fC'-^Ci _ 

{v, - Ei^^ + Ef^ + id){v, - 4"^ + 4 '^ + iS) 

7 -r( 0 )rr(l) rr(0)..(l) ^ 

_l_ 'ViM'v Mf'v fC'-^Ci 
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Here t'i/ = Vi — Vf, 4^^^ = 4^^ ~ 4^^ 4^^^ = 4^^ — 4^\ greater (lesser) projection of the molecular 

Green function, g^^^^ is greater (lesser) projection of Green function representing free electrons in nanoparticles 


9 k{T,T') = -iiTyCkir) cj'(r')), 


(S20) 
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and and G^^^^ are greater (lesser) projections of the mixed (molecule-nanoparticle) Green functions 


Gkx{T,T') = - i{TcCk{T) dlir')) 
Gxkir, t') = - i{Tc d^{T) c\{t')) 


(521) 

(522) 


In Eqs. (S16)-(S19) sum over s ( = {0,-,-!-}) corresponds to electronic (s = 0) and Stokes (s = —) and anti-Stokes 
(s = -|-) vibrational contributions to the total Raman signal with 


(S23) 



{Ny is average population of the vibrational mode). Note that Eq. (S19) contributes to vibrational Raman only. 
Its electronic part contributes to Rayleigh scattering [Ml |M] , and is disregarded in our consideration. Each of four 
expressions (S16)-(S19) contains the pure contribution of channel M in their first row (these terms represent the 
equivalent of model A, but without the effect of the e-v coupling (5) and the pure contribution of channel C in their 
second row. The rest of these expressions correspond to interference between these channels. It should be emphasized 
that each of these contributions has in turn a vibrational Raman and electronic Raman components, characterized 
by the corresponding change (vibrational or electronic) in system state. 

Numerical simulations presented in Fig. 5 of the paper show that this model is also capable of providing an 
asymmetric feature in the Stokes line. A closer examination of intra-molecular and charge-transfer contributions as 


well as their interference to vibrational and electronic Raman scattering (see Fig. S7) indicate that the Fano-like 
feature in the total signal is mostly due to charge transfer contribution to the vibrational Raman scattering. This 
contribution consists of multiple scattering paths from different electronic states in the contact and in this sense it is 
an interference feature. However, as in model A, the scattering lineshape asymmetry is dominated by the electronic 
Stokes sideband dressing the molecular vibrational line. 
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